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NUMERICAL  STUDY  OF  THERMAL  BLOOMING  WITH 
WIND  AND  TURBULENCE  AND  ITS  PHASE  COMPENSATION 

Zhang  Tianshu,  Lei  Guangyu,  Xie  Lijuan, 
and  Zheng  Shaotang 

Beijing  Institute  of  Applied  Physics 
and  Computational  Mathematics 

P.O.  Box  8009,  Beijing  100088  ^ 

*  ■  CT  In  this  paper  we  calculated  numerically  whole  beam  thermal  blooming  with 
.  lencc  and  their  ideal  phase  compensation  for  laser  atmosphere  propagation  usinj 
approximation  4D  code,  obtained  the  thermal  distortions  of  laser  beam  for 
<  id  analyzed  the  interation  of  thermal  blooming  with  turbulence.  In  addition  we  alsc 
le  selection  principle  of  time  step  and  space  step. 

ORDS  thermal  blooming,  turbulence,  phase  compensation. 


I .  Introduction 

There  are  various  favorable  and  unfavorable  effects  on  laser 
beam  transmission  due  to  atmospheric  conditions,  mainly  in 
variable  wind  force  and  variable  direction  relative  to  height  in 
the  atmosphere,  such  as  random  up-and-down  motion  of  wind  and 
turbulence  in  the  atmosphere.  Usually,  winds  are  favorable  to 
restrain  thermal  blooming;  however,  atmospheric  turbulence  mainly 
leads  to  the  effect  of  beam  expansion,  and  reduced  beam 
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coherence.  The  paper  presents  various  wind  models.  Turbulence 
is  simulated  with  multiphase  shield,  with  the  application  of  the 
Kolmogorov  spectrum.  In  section  II,  fundamental  equations  are 
given.  In  section  III,  the  selection  principle  of  time  and  space 
step  lengths  are  discussed.  In  section  IV,  the  results  of 
calculations  are  listed  along  with  a  brief  analysis.  In  section 
V,  several  conclusions  are  provided  in  the  stages  of  the  study. 

II.  Fundamental  Equations 

2.1.  Equation  of  the  Laser  Electrical  Field 

The  monochromatic  and  linearly  polarized  scalar  wave 
equation  is  obtained  in  the  paraxial  approximation 

-  +V? 
z 

In  the  equation,  v  ~  is  the  laser  wave  number.  In  the  paper, 
it  is  assumed  that  the  wave  numbers  for  the  main  laser  and  the 

beacon  light  are  the  same:  — - - 2-  .  The  laser 

propagation  direction  is  z .  6=1  corresponds  to  the  main  laser 

and  6=-l  corresponds  to  the  beacon  light.  The  change  6e  of  the 
dielectric  constant  in  the  atmosphere  is  due  to  thermal  blooming 
and  turbulence;  the  change 

Due  to  thermal  blooming,  the  variable  dielectric  constant 
SExb  is  proportional  to  the  variable  density  p^;  that  is, 

.46  p,  (x,  y,  z,  t) 
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Due  to  turbulence,  the  variable  dielectric  constant  SCxurb  is 
obtained  by  the  following  processing  of  the  equivalent  phase 
shield:  when  the  gap  between  the  turbulence  phase  shield  is  not 
large  (it  should  not  be  too  small),  otherwise  there  will  be  a 
statistical  correlation  between  the  two  phase  shields,  the 

f i^sction  effect  can  be  neglected.  Therefore,  from  Eq.  (1),  we 
obtain  the  field  equation  for  the  case  of  pure  turbulence, 
leaving  out  of  consideration  the  thermal  blooming  effect. 


zik 


dE 

dz 


+  A:  dc  tlirb 


Thus,  we  obtain 


ik 

^  ^  ^ ^  ^  ^  TURB 

In  the  equation,  AZp.,.  [p.s.=phase  shield]  is  the  distance 
between  the  phase  shields.  From  [1],  we  obtain 

i|  ^  d  Z 


I 


2n 


Az 


»» 


<I>  {K^,  )dKJK^ 


Here,  $  is  the  Kolmogorov  spectrum.  There  is 

In  the  equation,  the  turbulence  structure  constant  can  vary 

with  z.  x  ~  ,  Lq  is  the  external  dimension  of 

0 

turbulence.  a(K3t,Ky)  is  the  two-dimensional  complex  random 
function.  The  other  real  and  virtual  parts  are  subject  to  the 
Gaussian  distribution  with  1  as  variance,  that  is. 
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cos  2.  ' 


Re(a)=V  -2Inr, 

Im  ( (7 )  =  V  -  21n  Tj  S'"  2Jt 

In  the  equation,  .'(0,1)  is  the  pseudorandom  number.  To 

ensure  that  ,  is  a  real  number,  it  is  required  that  Re(a) 
and  Im(a)  satisfy,  respectively,  the  conditions  of  mirror 
symmetry  and  mirror-reflection  syitanetry. 

Eq.  (5)  explains  the  effect  on  turbulence  of  the  medium  in 
Azp  8  ;  this  is  equivalent  to  the  phase  shield  that  the  phase 
change  is  £  AZp,s.- 

In  the  case  of  pure  thermal  blooming,  the  field  equation  is 
similar  to  Eq.  (4),  with  the  only  difference  being  that 
is  replaced  with  len  .  The  solution  is  similar  to  Eq.  (5). 
Therefore,  only  the  phase  is  changed  in  the  pure  thermal  blooming 
that  neglects  the  diffraction  effect.  Only  when  considering  the 
diffraction  effect  (here,  the  term  V\E  is  included  in  the  field 
equation),  then  the  phase  change  induces  the  change  of 
oscillation  amplitude. 

2.2.  Linearized  Fluid-Mechanical  Equations 

By  linearizing  the  fluid-mechanical  equations,  the  following 

equations  of  mass,  momentum,  and  energy  are  obtained; 

+Po'^  •  =  0 

^  +Vp.  =  0 

dt 
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In  the  equations,  ...y,  and  px  are,  respectively,  the  varying 
quantities  density,  velocity,  and  pressure.  po  is  the 
atmospheric  temperature  without  turbulence,  c^  is  the  speed  of 
sound.  However,  is  the  ratio  between  the  specific  heats 

at  constant  pressure  and  at  constant  volume,  is  the 
absorption  coefficient.  The  equation  of  light  intensity  is 


££* 

The  total  derivative  is 


(c  is  the  speed  of  light) 


(11) 


In  the  paper,  only  the  transverse-direction  wind  is 
considered.  The  authors  used  the  four  wind  models  given  below: 

1.  x-direction  constant  wind:  = 

2.  x-direction  variable  wind:  it  is  assumed  that  there  is  a 

linear  relationship  in  ,  such  as 

Vx(z)=Vo+wz  (w  is  a  constant,  which  may  indicate  the 
rotational  angular  velocity  of  the  beam,  or  may  indicate  a 
constant  of  atmospheric  wind  field  varying  with  height). 

3.  Cyclostrophic  wind: 

We  assume  that  Vo(z)  is  a  constant  when  only  the  wind  direction 
changes  with  height  for  0(z)=|3zo. 

4 .  Random  wind 

r,  t)^V^(2)  +  syjx,y 
,  z.  t)=d  V,(x,y,  z,  t) 
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In  the  equations,  the  up-and-down  motion  (of  the  random  wind)  SV^ 
and  6Vy  are  calculated  with  simulation  processing  that  is  similar 
to  the  simulation  processing  of  SCturb* 

For  the  field  equation  ( 1 ) ,  the  authors  applied  two  methods 
for  its  solution.  One,  it  is  to  conduct  a  two-dimensional 
Fourier  transform  in  the  transverse  direction.  In  the 
longitudinal-direction  (propagation  direction),  the  phase  shield 
method  is  used  for  the  solution  of  two  levels  of  accuracy  [1]. 

In  the  second  method,  the  separation  of  components  is  used  to 
divide  Eq.  (1)  into  multiple-component  equations  [2].  With 
respect  to  transverse-direction  components,  one-dimensional 
Fourier  transform  is  conducted.  For  the  longitudinal-direction, 
the  analytical  solution  in  form  is  sought. 

The  isobaric  approximation  is  conducted  with  the  fluid- 
mechanical  equations.  Then  the  three  conservation  equations 


degenerate  into  an  equation  of  density  variation: 


,1,  I£l  +  y  l£i.  =  -  nl  a 
dx  ^  dy  c]  ‘ 


III.  Selection  of  Calculation  Region  and  Lattice 
3.1.  Selection  of  L  and  N 

To  calculate  for  the  functions  of  pseudo-transverse 
direction  wind,  and  to  reduce  the  boundary  effect,  the 
transverse-direction  period  L  should  be  as  large  as  possible, 
especially  during  large  variations  in  functions.  Since,  based  on 
the  periodic  boundary  conditions,  the  spectral  space  lattice 


interval  can  be  obtained  (we  make  ^  ).  It 

can  be  seen  that  only  when  larger  values  of  L  are  taken,  can  a 
small-spectral  lattice  be  obtained.  However,  the  maximum 

calculated  wave  number  .'?\=---[c  =  calculated]  (N  is  the 

...  / 

number  of  division  points).  We  can  see  that  if  the  high- 
frequency  portion  is  important  in  studying  the  problem, 
should  extend  to  the  high-frequency  band.  Thus,  when  L  is  large, 
N  is  also  large.  Selection  of  L  and  N  should  be  matched.  We 
take  L=6a  (a  is  the  beam  radius),  thus  (D=2a).  With 

respect  to  the  Kolmogorov  spectrum,  during  the  centimeter-level 
perturbation  wavelength,  the  value  of  the  spectral  function  has 
been  reduced  so  that  it  can  be  neglected.  Therefore,  when  D  is 
approximately  equal  to  Im,  N  should  be  greater  than  100  in  order 
to  let  cover  the  main  spectral  band.  Of  course,  with 

respect  to  the  problem  of  increasing  small  perturbations,  due  to 
the  large  gain  in  short-wave  perturbations,  this  is  the 
Kolmogorov  spectrum,  initially.  Then  very  quickly,  the  high- 
frequency  portion  will  increase,  thus  N  should  be  again 
increased.  However,  N  should  not  be  too  large  because  an 
overlarge  N  will  cause  large  calculation  errors  in  the  FFT 
method . 

3.2.  Selection  of  Longitudinal-Direction  Step  Length  Az 

When  using  the  multiphase  shield  method  at  accuracy  level  2, 
the  stability  requirement  is 

j  . 

5.8  X  10 -"Ip,  i* 
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when  A=l)im  and  under  the  isobaric  approximation,  there  is 

/■ 


Here,  MW/cm^  is  the  unit  for  Ip;  CGS  units  are  used  for  the  other 
quantities . 

With  respect  to  the  separation . of  components  method,  the 
longitudinal-direction  step  length  can  be  larger.  We  pick  At 
between  several  to  tens  of  meters . 

3.3.  Selection  of  Time  Step  Length  At 

With  respect  to  the  nonlinear  dynamics  problems,  it  is  best 
to  take  the  absolute  stability  format.  For  the  rapidly  varying 
stage.  At  should  be  as  small  as  possible.  Generally,  At  should 
be  between  0.1  and  0.01  of  the  processing  time.  We  select 
At=l/20(t„);  t„=a/Vo  is  the  wind  transition  time. 

IV.  Calculation  Results  and  Brief  Analysis 

In  the  calculations,  we  used  the  CGS  units.  The  simplified 
model  of  large  atmosphere  transmission  that  we  obtained  is  as 
follows:  the  homogeneous  atmosphere  at  sea  level  of  5km  thickness 
is  used  to  simulate  the  total  atmosphere  in  the  laser  path,  and 
we  assume  that  it  is  in  the  intermediate  turbulent  state.  The 
homogeneous  wind  velocity  is  500cm/s;  the  laser  wavelength  is 
selected  as  l^m;  and  the  laser  beam  radius  is  taken  as  50cm.  In 
the  paper,  the  beam  anomaly  compensation  adopts  the  ideal  phase 
conjugate  method. 


8 


4.1.  Beam  Rotation  (Equivalent  Shift  Shearing  Wind) 

Let  a  (rad/s)  be  equal  to  2.6xl0-^  5.2x10-3,  and  7.8x10-3  in 
three  states:  thermal  blooming,  thermal  blooming  with  turbulence, 
and  thermal  blooming  with  turbulence  and  phase  compensation,  in 
the  model  calculations.  Table  1  shows  the  calculation  results. 

For  few  models  in  Table  1,  we  compare  the  isobaric- 
approximation  results  in  this  research  and  the  non-isobaric- 
approximation  results  in  [3].  As  discovered  in  the  results  the 
difference  is  minor  at  three  times  the  wind  transition  time.  It 
is  apparent  that  isobaric  approximation  is  a  better 
approximation . 

From  Table  1,  when  Ip  is  large  (also,  the  value  of  the 
thermal  anomaly  number  is  large),  the  difference  between 
thermal  blooming  and  thermal  blooming  with  turbulence  is 
increased.  This  is  because  turbulence  will  induce  intensity 
fluctuations  61  (scintillation).  From  [4],  we  know  the  intensity 
fluctuation  scale  ‘ -  .  here  fj  is  the  focal  length  of 

thermal  blooming.  When  Ip  is  large,  fg  is  small,  thus,  I  is  also 
small.  However,  <  ^  ,  We  can  see  that  when  Ip  is  large, 

the  intensity  fluctuation  61  is  also  large.  From  Eq.  (1),  we 
know  that  large  61  values  will  also  induce  large  phase 
fluctuation  values  6$,  then  the  beam  expansion  effect  is 
intensified.  In  other  words,  the  larger  the  value  of  Ip,  the 
more  intensive  interaction  between  thermal  blooming  and 
turbulence.  Here  the  numerical  calculation  results  also  prove 
this  point. 
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«  .  Brightest  -point  Strehl  rttio  S  at  target  plane  for  three  conditions.  3r . ,  condition  A 
bloomin  gonly  B:  Thermal  blooming  with  turbulence,  and  C:  thermal  blooming  and  turbu 
with  ideal  cpmpensatlon  (the  laser  beam  intensity  /^Is  In  W/onS  the  slewing  rate  o)  In  ri. 


(o  N. 

Condi 

50 

100 

200 

400 

80(> 

A 

0.945 

0.89  (0.89) 

0.82(0.85] 

2.6  X  10-' 

B 

0.71 

0.55 

C 

0.91 

0.89 

0.79 

A 

0.98 

0.96 

0.92  (0.92) 

0.838(0.85  ] 

0.65 

).2x  10-* 

B 

0.76 

0.73 

0.68 

0.60 

0.51 

C 

0.96 

0.95 

0.93 

0.92 

0.82 

A 

0.984 

0.97 

0.935  (0.93) 

0.87(7) 

0.75 

1 

[0.95] 

10.93] 

bo 

X 

o 

B 

0.76 

0.74  j 

0.70 

0.54 

c 

0.952 

0,945 

0.93 

0.90 

0.88 

•  mber  in  [  ]  is  presented  by  [3]  ; 

umber  in  (  )is  obtained  at 


From  Table  1,  we  can  see  that  phase  compensation  degrades  in 
the  light-intensity  range  listed  in  the  table.  This  is  because 
this  range  lies  in  the  PCT  region.  From  [5],  we  know  that  the 
maximum  compensable  light  intensity  is 

In  the  equation  ~2r^  (^i 

transmission  system).  With  respect  to  our  model,  F=314, 

3W/cm^.  Therefore,  '  ,,  ';50W/cm^.  The  case  is  degraded 
with  more  compensation  when  exceeding  50W/cm^. 

We  know  from  the  data  listed  in  the  table,  the  beam  rotation 
angular  velocity  ta  varies  by  threefold,  but  the  variation  of  the 
Strehl  ratio  of  the  bright  spots  is  not  obvious.  This  is  so 


10 


because  the  propagation  distance  in  our  model  is  relatively 
short,  and  the  equivalent  shift  shearing  wind  is  not  intense. 
Furthermore,  we  can  see  that  the  turbulence  effect  is  also 
unrelated  to  the  beam  rotation  speed. 

From  column  3  of  Table  1  we  know  that  there  is  little 
difference  in  the  calculation  results  at  3t„  and  6t„.  Thus  we 
know  that  thermal  blooming  image  has  become  stabilized  at  3t„. 

4.2.  Cyclos tropic  Wind  (Only  Wind  Direction  Variation  is 
Considered  Here) 

In  the  calculation,  we  take  /3  (rad/cm)  equal  to  0.175xl0‘^ 
(50°  rotation  at  a  5km  height).  Ip  [W/cm2]=400.  Other  conditions 
have  been  stated  above.  With  respect  to  pure  thermal  blooming. 
Fig.  1  shows  the  calculation  results.  From  the  figure,  when  the 
image  rotates  by  a  certain  angle,  the  extent  of  thermal  anomaly 
does  not  change  because  the  thermal  anomaly  variation  number 

i,  -  .  With  the  wind  velocity  at  constant  and  only 

varying  the  wind  direction,  from  our  wind  model  we  know 

'  4.  j-  .  Here  Np  is 

unrelated  to  (i .  In  other  words ,  Nd  does  not  rely  on  the 
variation  of  velocity  or  direction. 

4.3.  Calculation  During  High  Thermal  Anomaly  Number 
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With  respect  to  CW  pulses,  under  the  homogeneous  medium 
conditions,  the  dimensionless  thermal  anomaly  number  is 

(here,  the  unit  of  Ip  is  W/cm^;  CGS 
units  are  used  for  the  other  quantities).  Here,  .  (a,  is 

the  scattering  coefficient).  For  our  model,  when  Ip=800W/cm^, 
Nd«2500.  Now  the  phase  difference  is  greater  for  adjacent 
lattice  points.  There  will  be  greater  error  in  calculations 
using  the  FFT  method.  Moreover,  there  are  perturbations  of 
varying  scales  at  this  time,  as  the  relative  importance  of  the 
high-frequency  band  of  the  spectrum  is  increased.  Therefore,  the 
number  of  lateral-direction  division  points  should  be  at  least 
greater  than  200.  In  our  calculations,  128  division  points  are 
applied.  So  there  will  be  relatively  larger  errors  for  this 
light  intensity.  (Refer  to  the  last  column  of  Table  1,  and 
Fig.  2.) 

V.  Conclusions 

Here,  the  listed  conclusions  are  the  research 
accomplishments  in  steps.  They  are  as  follows: 

5.1.  There  are  intensive  interactions  among  large  light 
intensity  (or  intensive  thermal  anomaly),  turbulence  and  thermal 
blooming.  This  results  of  scintillation  amplified  by  thermal 
blooming . 

5.2.  For  short  wavelengths  (iMm),  low  wind  velocity,  greater 
beam  radius  (in  the  vicinity  of  Im),  and  longer  propagation 
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distance  (longer  than  several  km)  can  compensate  for  low 
threshold  of  light  intensity  (below  lOOW/cm^) .  Exceeding  this 
threshold  value,  the  situation  will  be  degraded  with  more  phase 
compensation . 


\n  intensity  contour  at  3/^  ( /^*400W/cm^  Fig.  2  An  intensity  contour  at  3  u  whe  \ 

rotation  wind)  /, (inW/cm*)is  800 

5.3.  When  the  propagation  distance  is  not  long  (less  than 
several  km) ,  the  variation  in  beam  rotational  angular  velocity  is 
several-fold.  This  renders  little  effect  on  the  Strehl  ratio  of 
the  bright  points.  Furthermore,  the  effect  of  turbulence  is 
almost  unrelated  to  the  rotational  speed. 

5.4.  With  respect  to  a  CW  laser  beam,  the  thermal  blooming 
images  have  stabilized  at  threefold  the  wind  transition  time. 

5.5.  When  wind  velocity  is  constant,  and  wind  direction 
changes  only  with  height,  the  result  rendered  in  the  pure  thermal 
blooming  case  is  to  rotate  an  angle  for  the  entire  thermal 
blooming  image. 

5.6.  When  there  is  high  light  intensity  or  intensive  thermal 
anomaly,  it  is  difficult  to  calculate  thermal  blooming  by  using 
the  4-D  laser  transmission  program.  Now  there  is  physical 
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instability  and  also  numerical  instability;  therefore,  the 
horizontal  and  transverse-direction  lattices  should  become 
slender  and  match  to  each  other. 

The  authors  are  grateful  to  colleagues  Lin  Dewen  and  Xie 
Shumao  for  their  beneficial  comments  on  this  research. 

First  draft  was  received  on  June  1,  1993;  the  final  revised 
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CT  The  model  for  simulating  the  components  of  an  adaptive  optics  system  is  d< 
imerically  applied  to  the  laser  beam  propagation  through  turbulent  atmosphere 
phase  compensation.  Results  are  given  and  disccused  for  compensated  Streh 
with  atmospheric  coherent  length  anisoplanation  and  the  finite  bandwidth 
servo  system.  Also  the  variation  of  residual  phase  deviation  vs  djr^  is  given  ; 

-tRDS  adaptive  opticsi  numerical  simulation)  turbulence,  anisoplanation. 

I .  Introduction 

As  is  well  known,  there  are  many  successful  experimental  and 
applied  examples  of  phase  compensation  of  the  turbulence  effect 
on  laser  atmospheric  transmission.  As  revealed  in  recent 
research,  the  interaction  between  turbulence  and  thermal  blooming 
in  the  atmospheric  transmission  of  high-powered  lasers  will  lead 
to  instability  of  phase  compensation  [1,2].  In  fact,  how  many 
single  components  are  required  for  an  adaptive-optics  system  with 
a  finite-bandwidth  deformable  mirror?  With  what  kind  of  driver 
spacing,  system  response  bandwidth,  and  other  factors  can  one 
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attain  the  optimal  phase  compensation  effect  of  laser  atmospheric 
transmission?  This  series  of  problems  remains  to  be  answered. 
Thus,  a  numerical  analytical  model  and  a  computational  program 
are  established  for  numerical  simulation  of  the  adaptive-optics 
system;  this  is  important  in  studying  the  above-mentioned 
problems.  Furthermore,  this  is  an  indispensable  component  part 
in  the  four-dimensional  calculation  mode  of  laser  atmospheric 
transmission  with  compensation.  On  account  of  these  factors, 
many  scientists  made  outstanding  investigations  in  this  area 
[3,5].  The  paper  adopts  the  wavefront  detection  model  similar  to 
that  proposed  in  [4,5].  As  proposed  in  [3],  in  the  wavefront 
reconstruction  numerical  model,  the  wavefront  fitting  method  in 
the  practical  adaptive-optics  system  is  similar  to  the  method  in 
reference  [6].  On  this  basis,  the  computational  program  of 
numerical  simulation  was  compiled,  for  simulation  computation  on 
atmospheric  transmission  of  lasers  with  turbulence  under  phase 
compensation  of  the  adaptive-optics  system.  Under  identical 
compensation  conditions,  the  Strehl  ratio  at  the  far-field  of  the 
light  beam  and  the  residue  phase  variance  are  entirely  consistent 
with  the  simulation  results  in  reference  [3].  For  the  first 
time,  the  authors  conducted  computations  on  the  phase 
compensation  effect  by  finite  response  bandwidth  of  the  adaptive- 
optics  system  with  respect  to  different  transmission  aperture  and 
anisoplanatism.  The  computational  results  are  quite  consistent 
with  the  analytical  results.  This  explains  that  the  analytical 
results  of  the  simple  theoretical  results  is  reliable. 
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Fig*  1  A  schematic  diagram  of  adaptive  optics  system 


II.  Adaptive-Optics  System  and  Its  Numerical  Model 

Fig.  1  shows  the  principle  of  the  adaptive  optical  portion 
in  a  laser  transmission  system.  With  respect  to  the  anomalous 
wavefront  phase  of  real-time  measurement  of  the  beacon  light  with 
the  wavefront  detector,  the  detection  system  calculates  the 
corresponding  driving  quantity  of  the  deformable  mirror  driver 
and  the  corresponding  inclination  of  the  tilt  mirror,  which  is 
calculated  based  on  the  wavefront  detector  signal,  to  real-time 
drive  of  the  deformable  reflective  mirror  and  tilt  mirror,  thus 
producing  the  phase  wavefront  conjugated  with  the  beacon  light 
phase  to  apply  to  the  main  laser  in  order  to  compensate  for  the 
wavefront  anomaly  of  the  main  laser  caused  by  atmospheric 
transmission  during  turbulence  in  order  to  improve  light  beam 
quality. 

From  Fig.  1,  we  can  see  that  the  adaptive-optics  system 
includes  mainly  three  major  parts:  (1)  wavefront  detector;  (2) 
control  system;  and  (3)  tilt  mirror  and  deformable  reflective 
mirror  (which  is  the  wavefront  reconstruction).  In  the 
following,  the  mathematical  models  and  algorithms  of  these  three 
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portions  are  introduced. 


( 1 )  Wavefront  Detection 

Generally,  the  dynamic  Hartmann  method  is  used  for  wavefront 
detection.  The  light  spots  gravities  X(,„  and  y^n  of  the  beacon 
light  on  the  focal  plane  of  each  Hartmann  sublens  are  detected  by 
the  wavefront  sensor.  Thus,  the  corresponding  wavefront  tilt  of 
the  beacon  light  on  the  particular  Hartmann  sublens  can  be 
obtained. 

■  • 

In  the  equation,  F  is  the  focal  length  of  the  subaperture.  In 
numerical  simulation  the  light  spot  gravity  is  calculated  with 
the  following  equation  [4,5] 

Jim [  ]dr  jjj  \(f>Jdr 

n  indicates  the  n-th  sounding  sublens;  <|)b  is  the  field  function 
of  the  beacon  light;  k  is  the  wave  number  of  the  beacon  light. 

The  integration  region  is  the  region  of  the  sounding  subaperture. 

( 2 )  Control  System 

Mainly,  there  are  two  functions  for  the  control  system.  One 
is  to  use  the  wavefront  tilt  on  the  beacon  mark  thus  detected  to 
simulate  the  wavefront  of  the  beacon  light.  Then  the  driving 
quantity  of  the  driver  for  the  deformable  mirror  is  calculated. 
Secondly,  this  involves  carrying  out  the  hardware  driving  of  the 
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deformable  mirror.  Here  the  first  function  is  mainly  discussed. 
With  respect  to  the  second  function,  it  is  required  to  discuss 
with  the  information  control  theory  of  the  special  system.  In 
numerical  simulation,  the  simplest  method  is  to  introduce  the 
wavefront  reconstruction  error  brought  forth  with  the  control 
system  in  order  to  apply  to  the  conjugate  phase  for  processing. 

In  the  first  function,  the  wavefront  trial  fitting  can  adopt 
multiple  algorithms.  For  consistency  with  the  simulation 
algorithm  [6]  in  the  actual  adaptive  algorithm  system,  simulation 
is  carried  out  by  using  the  Zernike  polynomial.  In  this  method, 
first  the  wavefront  of  the  beacon  light  is  expanded  into  the 
Zernike  polynomial,  and  then  the  wavefront  tilts  on  the  various 
sublenses  are  detected  in  order  to  obtain  coefficients  of  various 
terms  of  the  expansion  formula. 

M 

Z  a,2,(,x,y) 

(-3 

M  is  the  number  of  expansion  terms  in  the  Zernike  polynomial. 

Then  the  mean  value  of  the  wavefront  inclination  on  the  plane  of 
the  n-th  sounding  sublens  is 

f  ^ 

,y)dxdyl A 

•  y 

Y.o,z,,(x  ,y)dxdyl  A 

The  integration  region  is  the  plane  of  the  n-th  sublens  and  its 
area  is  A;  the  subscripts  x  and  y  of  indicate  the  partial 
derivatives  of  x  and  y.  In  the  Zernike  expansion  equation. 
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summation  begins  with  the  third  term  (i=3).  This  indicates  that 
all  the  entire  tilt  items  of  the  beacon  light  wavefront  are 
separated  because  the  main  tilt  reconstruction  of  the  wavefront 
is  accomplished  by  the  tilt  mirror.  Similarly,  the  wavefront 
tilts  6x  and  0y  of  the  various  sublenses  thus  detected  are 
subtracted,  respectively,  from  the  total  average  tilts  63^  and  By 
of  the  wavefront.  Thus  we  obtain 

«-i  ii-i 


Let 


•  -I 


is  the  wavefront  tilt  variance.  Solve  for  'dalldoj  and  let 
•z  be  equal  to  zero,  and  then  we  obtain 

In  the  equation 


’  ,y)dxdylA 


,y)dxdylA 


[Zxin]  and  [Zyin]  are,  respectively,  the  transposed  matrixes  of 
[Zxn„]  and  [Zynm].  By  solving  the  linear  equation  (8),  we  obtain 
the  developed  coefficients  of  various  orders  of  Zernike  to  be 
substituted  into  Eq.  (3)  to  obtain  the  beacon  light  wavefront. 
This  wavefront  simulation  algorithm  is  called  the  mode  algorithm. 
We  also  conducted  experiments  on  other  wavefront  simulation 
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algorithms  (regional  method).  This  is  also  feasible  as  indicated 
in  the  results  of  numerical  calculations,  as  the  wavefront 
simulation  errors  are  basically  compatible  to  the  mode  simulation 
algorithm  of  the  Zernike  polynomial.  So  these  are  omitted  here. 
The  calculation  steps  for  the  driving  quality  of  the  deformable 
mirror  are  as  follows;  first,  the  wavefront  to  be  reconstructed 
is  expanded  into  the  linear  superposition  of  the  mirror  surface 
function  (which  is  also  called  the  influence  function)  of  the 
deformable  mirror  [3]. 


(10 


<-i 

Na  is  the  total  number  of  drivers  in  the  deformable  mirror,  d^ 
is  the  driving  quantity  of  the  i-th  driver.  Ii(x,y)  is  the 
mirror  surface  function  of  the  deformable  mirror  at  its  position. 
li(X/y)  can  use  the  measured  value  of  the  surface  type,  or  is 
described  in  the  form  of  functions .  When  describing  it  is  in  the 
form  of  functions,  generally  all  the  surface  type  functions  at 
the  driver  position  are  unified  into  a  function  form.  In  the 
actual  system,  of  course,  certain  errors  will  exist.  Here,  we 
apply  the  gaussian  function  to  describe  the  surface  of  the 
deformable  mirror. 


,  ,  ;  Xi  and  y^  are  the  position  coordinates 

of  the  i-th  driver;  d  is  the  mean  spacing  of  two  adjacent 
drivers;  p  is  the  coupling  coefficient.  In  other  words,  when  the 
driving  quantity  of  this  particular  driver  is  one  unit  and  when 
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the  drivers  at  the  adjacent  positions  do  not  act,  thus  it  is  the 
deformation  quantity  on  the  surface  of  the  deformable  mirror.  In 
the  following,  the  variance  <j\  between  the  reconstructed 
wavefront  and  the  beacon  light  wavefront  (simulated  wavefront)  is 
used  to  solve  for  the  reconstructed  beacon  light  <j)a  under  the 
least-squares  constraint: 

,  ^  •y)Vdxdy 


The  integration  region  is  the  entire  surface  of  the  deformable 
mirror.  Let  3  ,  we  obtain  [3] 

'  (r  (r  ,y)dxdy  =  ^^(p^in  ,y)Ij{x  ,y)dxdy 

0=1 ,2, -.AT,) 


When  written  in  the  matrical  form,  this  is 


In  the  equation. 


,(x  ,y)Ijix  ,y)dxdy 


l>,{x  ,y)lj(x  ,y)dxdy 


[lij]  and  [Rj]  are  called,  respectively,  the  influence  function 
matrix  and  the  simulation  matrix. 


(3)  Wavefront  reconstruction 

The  wavefront  reconstruction  of  the  tilt  mirror  in  the 
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numerical  model  is  obtained  by  integrating  the  mean  tilts  of 
Eg.  (10)  from  Eg.  (6).  The  wavefront  reconstruction  of  the 
deformable  mirror  is  obtained  by  substituting  into  Eg.  (10)  the 
driving  guantity  obtained  from  Eg.  (14).  In  the  actual  system, 
the  wavefront  reconstruction  principle  is  accomplished  by  the 
control  system  driving  hardware  (tilt  mirror  and  deformable 
reflective  mirror) .  Here,  naturally  a  certain  control  error 
exists.  As  revealed  by  analysis  in  reference  [3],  the  teinn  of 
the  largest  phase  compensation  error  in  the  adaptable  system  is 
from  the  wavefront  reconstruction  of  the  deformable  mirror. 
Generally,  this  is  several  times  to  tens  of  times  greater  than 
the  control  error.  We  can  see  that  the  effect  is  very  little  in 
the  results  of  numerical  simulation  by  neglecting  the  control 
error . 

III.  Results  of  Numerical  Simulation 

Based  on  the  above  discussion  on  model  and  calculation 
method,  a  niimerical  simulation  program  is  compiled,  to  conduct 
simulation  computations  to  compensate  for  the  turbulence  effect 
on  laser  atmospheric  transmission  in  a  37-component  adaptive- 
optics  system.  Under  the  turbulence  effect,  solving  of  the 
optical  wave  transmission  eguation  adopts  the  phase  shield 
method.  The  spectral  inverting  transmission  of  the  turbulence 
phase  shield  applying  the  folded-type  FFT  algorithm  as  the 
wavelength  of  the  beacon  light  is  the  same  as  the  main  laser 
wavelength  A=0 . 6328jjim) . 
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For  convenience  in  comparison,  the  same  parameters  in 
reference  [3]  are  selected  with  coupling  coefficient  p=0.15.  In 
the  numerical  calculation,  the  phase  distribution  of  the  beacon 
light  can  be  obtained  directly.  Thus,  from  Eqs.  (14)  to  (16), 
the  driving  quantity  of  the  driver  for  the  deformable  mirror  can 
be  solved.  In  other  words,  this  is  to  directly  apply  the  beacon 
light  phase  to  conduct  wavefront  reconstruction  calculations;  in 
other  words,  this  is  the  detection  and  trial-fitting  approach  to 
obtain  the  error  in  the  absence  of  a  wavefront.  The  wavefront 
reconstruction  residue  variance  and  the  Strehl  ratio  that  are 
obtained  are  consistent  with  the  results  of  numerical  simulation 
in  [3].  As  shown  in  solid  curves  in  Fig.  2a  and  b,  the  abscissa 
is  d/ro;  ro  is  the  transverse-direction  coherent  length  of  the 
atmosphere.  The  data  point  +  and  A (diamond-shape)  indicate, 
respectively,  the  numerical  simulation  results  when  the  emission 
apertures  D  are  1.67m  and  0.25.  Moreover,  consideration  is  given 
to  the  wavefront  detection  and  wavefront  trial-fitting  errors. 

In  other  words,  the  wavefront  inclination  is  calculated  with  the 
beacon  light  field  by  using  Eqs.  (1)  and  (2).  Then,  from  Eqs. 
(3),  (8),  and  (9),  the  wavefront  reconstruction  calculation  is 

conducted  after  wavefront  trial  fitting.  From  Fig.  2,  we  can  see 
the  following:  (1)  by  leaving  out  of  consideration  on  the 
wavefront  detection  and  trial  fitting  error,  this  is  obviously  an 
overestimate  of  the  effect  of  phase  compensation,  especially  when 
d/ro  is  greater.  When  d/ro=3.0,  the  wavefront  detection  error 
and  wavefront  trial-fitting  error  are  only  smaller  by  about  one- 
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ulalion  ;)•=  1.67m  ;  A  simulalton  0.25m  +  simulation  £»=  1. 67m  ;  O  simulation  O-=0.: 
considering  error  of  wave -front  reconstruction  (b>wave-  front  reconstruction  error  (numerical  : 

Fig. 2  Strehl  ratio  (a) and  residual  phase  variance  vs  d/r, 

half  of  the  wavefront  reconstruction  error.  (2)  Similarly,  under 
the  condition  of  the  same  d/ro,  but  different  emission  aperture, 
although  there  is  no  obvious  difference  in  the  residue  phase 
variance  (as  shown  in  Fig.  2b),  the  phase  compensation  effect  is 
different,  especially  so  in  the  case  of  more  intensive  turbulence 
(larger  d/ro).  Possibly,  the  main  reason  is  as  follows:  when  D 
is  smaller,  under  the  condition  of  intensive  turbulence, 
fluctuation  of  the  oscillation  amplitude  will  obviously  reduce 
the  compensation  efficiency  of  the  AO  system.  However,  when  D  is 
larger  but  d/ro  remains  unchanged,  turbulence  is  weaker.  Thus, 
together  with  the  aperture— smoothing  effect,  fluctuation  of 
oscillation  amplitude  will  be  reduced.  Therefore,  the  effect  of 
AO  compensation  is  better. 

Fig.  3  shows  the  effect  of  anisoplanatism  with  regard  to 
phase  compensation  efficiency  of  the  AO  system.  In  the  figure,  6 
twenty- two  spaces  is  the  isoplanatic  angle  [8];  0  is  the  included 
angle  between  the  beacon  light  and  the  main  laser;  the  solid 
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curve  indicates  the  results  of  the  theoretical  analysis  [8];  and 
the  longitudinal  coordinates  indicate  the  relative  Strehl  ratio . 
This  is  the  result  after  normalizing  the  Strehl  ratio  when  0=0 . 
Thus,  it  is  convenient  to  compare  with  the  theoretical  result 
[8].  Since  no  consideration  was  given  in  the  theoretical 
analysis  to  the  effect  of  the  number  of  limited  components  and 
fluctuation  of  oscillation  amplitude,  this  the  entire  phase 
conjugate  compensation.  From  the  figure,  the  numerical 
simulation  result  matches  very  closely  with  the  theoretical 
analytical  result.  When  e/0o=l.O,  the  relative  Strehl  ratio  is 
approximately  0.78.  In  other  words,  in  these  critical 
isoplanatic  conditions,  the  compensation  efficiency  of  an  actual 
system  is  lower  by  22%  than  the  compensation  efficiency  of  the 
entire  phase  conjugation  (0=0)  at  the  instantaneous  response.  In 
the  condition  of  different  d/r®,  the  relationship  between  the 
relative  Strehl  ratio  and  0/0o  is  basically  consistent. 

To  conduct  a  simulation  analysis  on  the  effect  of 
compensation  efficiency  by  time  response  bandwidth  in  an 
adaptive-optics  system,  we  should  consider  the  turbulence  phase 
shield  varying  with  time.  We  know  that  under  the  approximate 
conditions  established  by  Taylor's  assumption,  the  variance  of 
turbulence  with  time  can  be  converted  into  variation  of 
translation  of  wind  speed  with  space.  Therefore,  we  will  apply 
the  translation  of  wind  speed  with  the  turbulence  phase  shield  to 
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Fig.  3  Relative  Strehl  ratio  0/0,  Fig.  4  Strehl  ratio  vj  servo  bandwidth  at  i/ 

m  3  S  4  jKl%^%«KXd-Strehl 

-theory:  +  simulation  dlr,=  l.O  — Theory  ;□  simulation  D=l.0m 

O  simulation  ■*■  simulation  D”0.48m 

simulate  variation  with  time  [4].  Besides,  it  is  still  generally 
true  that  the  response  bandwidth  of  the  adaptive-optics  system 
can  be  simplified  into  a  rectangular  function.  In  other  words, 
this  corresponds  to  the  time  lag  At=l/fo.  fo  is  the  response 
interruption  frequency  of  the  adaptive-optics  system.  Fig.  4 
shows  the  comparison  of  the  results  between  numerical  simulation 
and  theoretical  analysis  [8].  In  the  figure,  the  relative  Strehl 
ratio  shown  is  the  Strehl  ratio  under  the  instantaneous  response 
compensation  conditions .  The  Strehl  ratio  is  then  normalized  as 
a  vector.  The  purpose  is  to  eliminate  the  effect  on  the  Strehl 
ratio  due  to  residue  phase  difference  in  the  actual  system 
because  of  oscillation  amplitude  fluctuations  and  limited 
components.  This  effect  has  been  described  in  Fig.  2.  Thus,  the 
effect  on  the  Strehl  ratio  due  to  the  system  response  bandwidth 
is  revealed.  In  the  figure,  nug  indicates  the  mean  wind  speed; 
the  value  of  d/rg  is  1.0.  The  solid  curve  in  the  figure 
indicates  the  result  of  theoretical  analysis  (in  the  case  of  pure 
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conjugate  phase  compensation).  We  can  see  that  the  result  of 
numerical  simulation  after  normalization  is  quite  consistent  with 
the  theoretical  analytical  results.  In  other  words,  it  is 
rational  to  use  the  wind  speed  translation  model  to  describe  the 
time  correlation  of  turbulence.  Besides,  if  the  space 
correlation  scale  l=2DELTAx  is  used  for  the  turbulence,  the  scale 
is  the  reciprocal  number  of  the  Nyquist  sampling  frequency,  and 
Ax  is  the  transverse-direction  calculated  sampling  width.  Under 
Taylor's  assumptions,  the  characteristic  frequency  (of 
turbulence)  is  fx=  ^  .  Then  when  fo=fT/  the  relative  Strehl 

ratio  is  approximately  0.88,  which  is  wholly  consistent  with  the 
value  from  theoretical  analysis  [8].  In  other  words,  by  using 
the  space  relationship  of  turbulence  to  describe  its 
characteristics,  this  provides  a  convenient  and  feasible  means  in 
studying  the  effect  on  compensation  efficiency  due  to  the 
response  bandwidth  of  the  adaptive-optics  system. 

IV.  Conclusions 

In  this  paper,  the  numerical  model  and  algorithm  of  an 
actual  adaptive-optics  system  are  described.  Under  the  same 
calculation  conditions,  the  numerical  calculation  results  fully 
match  the  results  from  numerical  simulation  in  [3].  The  effects 
on  compensation  efficiency  by  different  emission  apertures  D, 
atmospheric  coherent  length  Tq,  and  anisoplanatism,  the  system 
limited  response  bandwidth  is  very  consistent  with  the  results 
predicted  from  theoretical  analysis.  With  regard  to  other 
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factors  of  influence  in  the  adaptive-optics  system,  such  as 
coupling  coefficient,  mirror  surface  response  function,  and 
variations  in  the  inhomogeneous  atmospheric  conditions  in  the 
actual  case,  the  analysis  of  compensation  efficiency  in  the 
adaptive-optics  system  is  a  problem  requiring  further  and  deeper 
research.  However,  the  nvimerical  calculation  with  compensation 
experiments  in  the  actual  system  will  be  an  important  means  of 
further  improving  the  simulation  model.  Moreover,  in  the 
theoretical  models  [ 6 ]  in  the  paper . and  the  other  models 
compared,  it  is  simpler  to  describe  the  adaptive-optics  transfer 
function.  Therefore,  eventually  the  experimental  results  should 
be  relied  on  when  examining  the  numerical  model.  This  is  the 
next  step  in  research. 

The  first  draft  of  the  paper  was  received  on  April  20,  1993; 
the  final  revised  draft  was  received  for  publication  on  July  16, 
1993. 
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lACT  The  optical  performance  of  trichromatic  beam  splitters  developed  for  DF  last. 

as  shared  -aperture  components  is  presented  .  In  this  paper  .  we  have  discus 
ite  of  splitter  optical  elements  meeting  spliting  performance  .the  selection  of  coating 
’".ultilayer  film  design  and  techniques  .and  the  performance  results  measured  .The  exp 
Jits  of  R>99.6%  at  3.8/im  wavelength  .the  average  transmittance  of  more  than  0.7: 
avelcngth  8  ~  14^m  and  the  average  transmittance  of  more  than  0.50  at  the  CGD  rt 
for  the  trichromatic  splitting  film  with  uncoated  back  surface  of  ZnSe  substrate 
.  The  environmental  stability  of  the  samples  are  examined  . 

^ORDS  DF  laser .  share  -aperture  components,  trichromatic  coating,  CVD  ZnSc  sub. 

I .  Introduction 

High-performance  trichroic  beam  splitters  are  used  in 
shared-aperture  infrared  DF  (deuterium  fluoride)  chemical  lasers. 
The  beam  splitter  requires  low-absorptance  film-coating  material. 
Such  shared-aperture  device  has  a  total-reflection  DF  laser 
wavelength  band.  Moreover,  the  splitter  has  dual  functions  of 
easily  transmitting  CCD  spectral  response  and  a  long-wave 
infrared  spectrum  between  8  and  14)im.  The  shared -aperture  system 
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has  the  advantages  of  higher  resolving  power  and  signal-to-noise 
ratio.  Moreover,  the  shared  optical  path  in  the  system 
simplifies  optical  collimation  and  tracking.  However,  it  is  very 
difficult,  experimentally,  to  share  the  multichroic  film  of  the 
shared-aperture  optical  component.  J.  E.  Rudisill  et  al. 
fabricated  a  dichroic  beam  splitter  with  reflection  at  3.8|im  and 
transmittance  between  8  and  144m  [1].  To  fabricate  the  trichroic 
beam  splitter,  it  is  required  co  deposit  a  multilayer  film  medium 
on  a  highly  transparent  substrate  from  the  visible-light  waveband 
to  the  infrared  waveband  (14|im).  To  upgrade  the  optical 
performance  and  the  film  layer  system,  the  absorption  inside  the 
substrate  and  at  the  surface  should  be  as  small  as  possible. 

The  article  discusses  the  design  and  preparation  of  the 
shared-aperture  trichroic  beam  splitter.  In  experiments  on  the 
trichroic  beam  splitter  film,  there  is  reflectivity  of  greater 
than  99.6%  at  the  3.84m  DF  output  band.  In  the  CCD  optical 
region  between  0.64  and  0.924m,  the  mean  transmittance  is 
approximately  50%;  between  8  and  144m  of  infrared,  the  mean 
transmittance  is  greater  than  75%. 

II.  Selection  of  Substrate  and  Film-Coating  Materials 

As  a  substrate  material  for  the  shared-aperture  device  in 
the  DF  laser  system,  it  was  required  to  have  high  transparency  in 
the  visible-light  region  and  the  long-wave  infrared  region,  in 
addition  to  good  mechanical  strength.  There  are  very  few 
materials  capable  of  transmitting  in  the  visible  light  region  and 
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also  transmitting  in  the  window  (8  to  14|jim)  infrared  region.  The 
most  commonly  used  materials  are  ZnSe,  ZnS,  BaFj,  alkali  metal 
halides,  and  diamond.  Although  alkali  metal  halides  have  the 
optical  properties  of  high  transparency  from  the  ultraviolet  to 
the  long -wave  infrared  region,  still  their  mechanical  strength 
and  environmental  stability  are  very  poor.  Diamond  is  expensive 
and  large  diamonds  are  even  rarer.  Transparency  in  the  visible 
light  region  for  large  CVD  ZnS  is  low;  its  mean  transmittance  is 
lower  than  40%  between  0.5  and  l.Ojim.  Moreover,  in  the  visible 
light  region,  the  transmittance  curve  varies  with  inclination 
[2].  From  0.6  to  17}im,  CVD  ZnSe  is  highly  transparent.  BaF2  is 
transparent  up  to  124m,  then  transparency  begins  to  decrease.  At 
14|im,  transparency  is  chancy. 

For  chemical  vapor  phase  deposited  large  ZnSe  pieces. 

Table  1  shows  their  mechanical  properties.  Fig.  1  presents  a 
curve  showing  the  CVD  Raytran  ZnSe  optical  spectrum  as  plotted 
experimentally  at  Raytheon  Corporation  in  the  United  States.  For 
^^^5®  pieces  of  ZnSe  material ,  its  absorptance  is  very  low; 
between  3.8  and  5.25^m,  its  absorption  coefficient  is 
4xl0"^cm"^ . 

For  BaF2  and  ZnSe,  both  materials  can  be  used  as  the  window 
material  of  the  shared-aperture  window  optical  device.  In  the 
authors'  research,  ZnSe  was  used  as  the  substrate.  In  the 
shared-aperture  spectroscopic  device,  the  substrate  should  be 
highly  transparent  in  the  visible  light  region  and  in  the 
infrared  region.  In  addition,  its  film-coating  material  should 
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Table  I  Typical  properties  of  CVD  ZnSe’"’ 


g.'  cm’) 

5.27 

index  of  refraction  inhomogencity 

.r  (^m) 

70 

(ppm.  max ) 

trength  (psi, 

temperature  cocfricjcnt 

l  lobding ) . 

7500 

(dn;  dT  *C  )at  I0.6^im 

modulus  (psi) 

9.75  X  10* 

extinction  cocfncicni 

;  ratio 

0.28 

632.8nm  tem'*) 

(Knoop  50gm ) 

100 

(g  1.06pm  (cm") 

/  (n-  cm  ) 

-  10’^ 

bulk  absorption  coefficient 

•:x  pa  ns  ion  /  t) 

@  2.77pm  (cm") 

:  20~170t: 

7.57  X  10'* 

@  3.8pm  (cm") 

•  i:orductivity  (251C,cgs) 

0.043 

@  5.25pm  (cm") 

cat  (cal/gm  T  ) 

0.081 

10.6pm  (cm'') 

transmission 

>70% 

pulse  damage  threshold. 

.;ion  limits 

0.5  ~  22/im 

peak  intensity  (GW  cm’) 

.  index 

@  180  ns 

3  ^m 

2.417 --  2.385 

(§  2.7pm 

/Lim 

2.403 

@  10. 6pm 

also  be  highly  transparent  in  both  regions .  At  the  wavelength 
3.85jiin  of  the  laser  device,  its  absorptance  is  very  low.  Fig.  2 
shows  the  absorptance  at  3.8(jini  for  some  film  layer.  Thus,  ThF/^, 
PbF2,  and  YbFj  can  be  selected  as  the  film-coating  materials. 
Since  NaCl,  KCl,  and  BaFz  are  easily  softened  by  absorbing 
moisture  and  since  YbFs  film  is  very  highly  stressed,  these 
materials  are  not  readily  mutually  compatible.  Finally,  the 
highly  refractive  film-coating  materials  ZnSe  and  ZnS  were 
chosen. 

III.  Design  of  Multilayer  Film 
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n  0.6  0.8  1.0  2.0  8.0  10 

Wavele..^ 

Fig.  I  Typical  transmittance  of  C  (20mm  thick) 

Any  periodic  multilayer  film  has  a  reflection  band  and  a 
transmission  band.  The  necessary  condition  for  a  reflection  band 
is 

In  other  words,  the  summation  of  the  optical  thicknesses  of 
various  film  layers  in  the  fundamental  period  are  1/2  of  the 
integer  times  the  center  wavelength  In  the  equation,  n^di  is 

the  optical  thickness  of  the  i-th  layer  in  the  fundamental 
period;  q  is  the  integer.  Aq  is  the  center  wavelength  possibly 
appearing  in  the  reflection  band.  However,  the  above-mentioned 
conditions  are  insufficient.  For  certain  wavelengths,  although 
the  summation  of  film  thicknesses  in  the  fundamental  period  is  a 
multiple  of  1/2,  yet  only  when  the  thickness  of  various  layers 
are  an  integer  times  1/2  the  wavelength  Aq,  can  the  various  film 
layers  seem  to  disappear,  as  is  also  the  case  for  the  fundamental 
period  and  the  entire  film  system  with  respect  to  the 
transmission  band  corresponding  to  these  wavelengths.  Therefore, 
based  on  Eq.  (1),  the  center  wavelength  of  the  reflection  band 
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can  be  determined,  and  the  transmission  band  can  also  be 
determined.  For  the  shared-aperture  trichroic  spectroscopic 
film,  it  is  required  that  the  reflection  band  be  at  the  center 
wavelength  Ao=3.8jim.  Moreover,  it  is  required  to  have 
transmission  in  the  spectral  response  region  (6.0  to  6.90nm)  for 
the  CCD  camera.  In  other  words,  this  is  required  that  the  center 
range  of  the  wavelength  Ao/5=0.76|im  be  in  the  transmission  band. 
With  respect  to  the  standard  quarter-film  stack,  this  wavelength 
region  is  the  reflection  band,  not  rip  to  the  requirements  of  the 
spectroscopic  properties.  If  the  standard  quarter-film  stack  is 
not  used,  the  nonstandard  film  system  can  result  in  a 
transmission  band  between  the  0.76jim  center  region  and  the 
infrared  region  between  8  and  14|j,m,  while  the  3.8|im  wavelength 
region  is  the  reflection  band.  For  simplicity  in  preparation, 
and  as  few  as  possible  film-coating  materials  to  be  used,  two 
film-coating  materials  are  used  to  design  the  film  system  in 
achieving  the  spectroscopic  properties.  Following  this  design, 
film  systems  of  the  four  following  periods  can  meet  the 
requirements  of  spectroscopic  properties:  4  32s.  H  represents  a 
ZnSe  film  at  one-quarter  of  refractance  at  3.8|am  wavelength;  L  is 
the  PbF2  or  ThF/;  film  at  low  one-quarter  refractance. 

The  sximmation  of  the  optical  thicknesses  of  the  various  film 

layers  in  the  fundamental  period  is 

'  '  *1  1 

With  the  reflection  band  L  appearing,  the  center  wavelength 
satisfies  the  condition 
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With  respect  to  g=5,  10,  that  is,  wavelengths  of  Xq/S, 

Aq/IO,  ...,  each  film  layer  seems  to  be  nonexistent.  Therefore, 
the  multiple  film  layers  in  the  fundamental  period  and  composed 
of  the  fundamental  period  seem  to  be  nonexistent.  Therefore, 
these  wavelengths  corresponds  to  the  transmission  band. 

Obviously 

g=l,  2,  3,  4,  6,  7,  ...,  reflection  band 

g=5,  10,  ...,  transmission  band 
In  other  words,  in  the  region  Ao=3.8)im  wavelength  is  the 
reflection  band;  and  the  region  Ao/5=0.76nm  corresponds  to  the 
transmission  band.  For  the  absorption  rate  in  the  DF  laser 
waveband  to  be  lower  than  0.05%,  the  film-coating  material  was 
selected  so  that  it  has  low  absorptance  at  the  3.8|im  wavelength, 
thus  reducing  the  thickness  of  the  entire  film  system.  For 
reflectance  in  the  3.84m  wavelength  band  to  be  at  the  maximvim  and 
transmittance  to  be  reduced  to  less  than  0.1%,  with  the  minimum 
number  of  film  layers,  it  is  required  that  both  film  coating 

materials  have  a  refractance  ratio  as  high  as  possible.  With 

consideration  given  to  the  laser  damage  threshold  value,  the 
outermost  layer  should  be  the  low-absorption  layer;  its  electric 
field  intensity  should  be  low.  Therefore,  two  kinds  of  film 

structures  7  r  /  and  -  ,  can  be  used  as  the  fundamental 

period  of  the  multilayer  film  for  the  shared -aperture  device.  By 
using  Ao=0.76|Jim  as  the  reference  wavelength  for  the  design. 
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1.2  3. S/iin  absorption  in  single-layer  1 

on  CaF,  substrates"' 

while  the  fundamental  period  becomes  4L6H  and  6H4L,  the  12-pair 
multilayer  film  structure  Sub/ ( 4L6H) ^^/Air  (composed  of  PbF2  and 
ZnSe  on  a  ZnSe  substrate)  should  have  its  theoretical  spectrum 
response  curves  as  shown  in  Fig.  3;  the  incidence  angle  of  the 
light  beam  is  22.5°. 

Fig.  3a  indicates  the  transmittance  of  the  wavelength 
response  in  the  DF  laser.  Fig.  3b  and  3c  show  the  CCD  spectral 
responses  and  the  8  to  14|im  transmittance  curves.  With  respect 
to  the  wavelength  between  8  and  14}i.m,  we  do  not  expect  that  the 
theoretical  calculations  are  accurate  because  the  optical 
constants  of  the  film  coating  materials  in  this  spectral  region 
are  uncertain. 

In  Fig.  3b  and  3c,  the  spectral  response  does  not  include 
the  effect  of  reflectance  at  the  back  of  the  ZnSe  substrate 
(theoretically,  over  17%)  of  the  beam  splitter. 

At  the  back  of  the  beam  splitter,  we  should  deposit  a 
reduced-reflectance  film  with  two  widely  separated  wavebands 
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Fig,  3  Theoretical  spectral  response  for 
(PbF,/ZnSc)‘*  on  ZnSc  at  22.5  * 

(a) 3.8^m  region  transmittance 

(b) CCD  response  region  transmittance 
(c  )8  —  I4^m  response  transmittance 


between  the  visible  light  region  and  the  infrared  region  between 
8  and  14inn/  in  order  to  eliminate  the  reflectivity  wear  onback. 
Table  2  shows  a  trichroic  spectroscopic  film  Sub/ ( 4L6H)^^/Air  in 
a  single-surface  ZnSe  substrate;  at  the  back,  there  are  the 
expected  properties  of  the  reduced-reflectance  film.  As  shown  in 
Table  2,  the  required  properties  can  be  satisfied  by  combining 
these  two  kinds  of  films. 
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Table  2  Theoretical  coating  designs  for  22.5^  incident  on  ZnSe 
for  a  3.8/jm  reflecting  trichroie  beam  splitter 


DF  laser  respoase 

Hector 

optical 

3.6  —  4.2^m  average 

mating 

thickness 

reflectance  transmiiiance 

absorptana 

ssign 

at  3.8/ini 

(%)  (•/.) 

(%) 

•<F,/ZnSef 

(0.20/0.30) 

99.92  0.03 

0.05 

long  waveiengtli  band  and  visible  CCD  response 

AR 

optical 

8  ~  14;^  average 

0.64  — 0.92;jim  ave 

coating 

thickness 

transmittance*  {%) 

transmittance* 

design 

at  0.76^m 

goal  predicted 

goal  predicte 

iS/PbF,/BaF,) 

(0  25/3.0/0.25) 

80  88 

SO  88 

*  Total  includes  both  surfaces 


IV.  Experiments 

To  prepare  a  high-quality,  low-absorptance  surface,  carving 
and  polishing  were  used  to  grind  and  polish  the  ODSOmm  CVD  ZnSe 
so  that  a  good  (0.1  diaphragm)  and  0.5'  parallelism  can  be 
obtained  at  the  substrate  surface.  The  post-polishing  substrate 
was  washed  with  a  mixture  of  ethyl  alcohol  and  ethyl  ether;  after 
drying  with  nitrogen,  the  substrate  was  immediately  placed  into  a 
film-coating  chamber.  Prior  to  film  coating,  an  Ar*  ion  beam  was 
used  to  bombard  the  substrate  for  20min  to  further  remove  dirt 
from  the  surface.  The  ion  beam  energy  was  300ev  and  the  beam 
current  density  was  20mA/ cm^. 

The  thin-film  deposition  was  done  in  an  AQ700  film-coating 
machine;  gases  were  evacuated  from  the  machine  with  a 
turbomolecular  pump  so  that  the  vacuum  was  free  of  contamination. 


40 


A  spherical-surface  rotating  jig  was  used  to  improve  film  layer 
homogeneity .  The  degree  of  vacuum  with  the  instrument  was 
lO'^^a.  There  was  a  hot-cathode  ion  source  installed  in  the  film¬ 
coating  machine  for  ion-beam  supplementary  deposition.  Both  PbF2 
and  ZnSe  were  deposited  with  the  electron  beeim.  The  PFj 
deposition  gas  pressure  was  6.2xlO'^Pa  and  the  ZnSe  gas 
deposition  pressure  was  3.4xlO'^Pa.  The  substrate  temperature 
was  160°C.  The  film  thickness  was  controlled  with  an 
optoelectronic  limit  value  method,  with  a  control  wavelength  of 
0.764m.  A  quarter-layer  of  ZnSe  film  at  0.764m  was  coated  in 
advance  on  a  K9  glass  control  plate;  the  method  of  multiple 
limiting  values  was  used  to  control  the  spectroscopic  film. 

V.  Results  of  Experiments 

On  the  front  surface  of  the  ZnSe  substrate,  a  trichroic 
spectroscopic  film  of  reflection  at  3.84m  (as  mentioned  in  Table 
2)  was  deposited.  By  using  a  Perkin-Elmer  350  infrared 
spectrophotometer,  the  infrared  transmission  spectrum  of  the 
coated  specimen  was  measured.  By  measuring  the  central  position 
of  the  reflection  band,  the  film  layer  thickness  was  examined. 

By  adjusting  the  control  wavelength,  the  center  of  the  reflection 
band  was  placed  at  3.84m.  Since  the  reflection  band  width  is 
difficult  to  measure  the  position  of  maximum  reflectance  in  the 
reflection  band,  the  authors  measured  two  wavelengths  when  the 
transmittance  rose  to  1%.  The  mean  value  of  these  two 
wavelengths  was  defined  as  the  center  wavelength. 
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Fig.  4  Measured  spectral  performance  of  ZnSe  3.8;cTn  reflecting  trichroic  beam  splitter 
a .  3.8;xm  band  and  8  ~  Uyim  long  wavelength  performance 
b.  (1  )0.5  visible  spectral  performance  ;  (2 ) transmittance  of  uncoated  ZnSe  substra 


Film  thickness  is  controlled  with  multiple  limiting  values 
by  using  transmitted  light.  The  controlled  wavelength  was 
determined  by  a  grating  monochromator.  After  multiple 
experiments,  consideration  was  given  to  chromatic  dispersion  of 
the  film-coating  material,  as  well  as  the  difference  in  the 
deposition  thickness  between  the  control  substrate  and  the 
specimen;  finally,  the  film  thickness  was  controlled  at  0.76)im 
wavelength.  Fig.  4  indicates  the  spectral  experimental 
properties  of  the  coated  specimens  in  the  infrared  zone  and  the 
visible  light  zone.  As  shown  in  the  figure,  under  the  condition 
that  the  reduced-reflectance  film  is  not  coated  on  the  back  of 
the  ZnSe  substrate,  and  the  high-reflectance  band  is  in  the  range 
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between  3.4  and  4.2|im  wavelength,  with  mean  transmittance  less 
than  0.4%,  the  measured  reflectance  at  3. Sum  is  greater  than 
99.6%;  and  the  mean  transmittance  in  the  spectral  region  between 
8  and  14^m  is  greater  than  75%;  and  the  mean  transmittance  in  the 
visible-light  spectral  region  between  0.64  and  0.92^m  is  close  to 
50%.  If  a  reduced-reflectance  film  with  the  widely  separated 
dual  wavebands  is  coated  on  the  back  of  ZnSe  as  mentioned  in 
Table  2,  the  transmittance  properties  in  the  visible-light 
spectra  and  the  spectral  region  between  6  and  14|jim  will  be 
further  improved. 

While  making  the  optical  property  measurements  on  the  film 
layers,  the  authors  evaluated  the  environmental  stability  and  the 
adherence  property  of  the  film  layers.  The  stability  experiments 
were  conducted  as  per  International  Standard  GB  1317-77.  The 
film  layer  experiments  involved  moisture  in  the  air,  low 
temperature,  and  saline  water.  After  these  experiments,  the 
optical  properties  of  the  film  layers  remained  unchanged.  With 
adherence  force  examination  of  the  film  layers  by  using  adhesive 
tape,  the  adhesive  tape  was  unable  to  separate  the  film  layers. 

The  authors  conducted  two  coating  methods  on  the  film 
deposition.  One  involved  general  thermal  evaporation.  In  the 
other  method,  while  thermal  evaporation  took  place,  Ar^  ion-beam 
bombardment  was  conducted;  the  bombardment  beam  current  was  60mA. 
It  was  discovered  that  under  a  spotlight  with  high  transmission 
in  the  visible-light  region  and  high  thermal  evaporation  in  the 
trichroic  spectroscopic  film  with  ion-beam  supplementary 
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deposition,  the  optical  scattering  was  obviously  less  than  that 
for  the  thermally-evaporated  film.  This  is  possibly  explained  by 
the  fact  that  the  bombardment  film  was  denser  than  the  evaporated 
film,  and  the  coarseness  was  less  than  that  of  the  evaporated 
film. 

The  first  draft  of  the  article  was  received  on  June  1,  1993; 
the  final  revised  draft  was  received  for  publication  on  July  30, 
1993. 
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